Impurity bands and sequential resonant tunneling in the presence 

of terahertz fields 



Andreas Wacker and Antti-Pekka Jauho 
Mikroelektronik Centret, Danmarks Tekniske Universitet, DK-2800 Lynghy, Denmark 



Stefan Zeuner and S. James Allen 
Center for Terahertz Science and Technology, University of California at Santa Barbara, Santa 

Barbara, California 93106 
(email: wacker@mic.dtu.dk) 



Abstract 

A theoretical and experimental study of transport in a low doped multiple 
quantum well structure shows that impurity bands are essential in understand- 
ing the electronic transport both with and without terahertz irradiation. A 
full self-contained model, which involves only the nominal parameters of the 
sample, agrees quantitatively with a wide range of experiments. 
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Perpendicular charge transport in biased superlattices displays a rich variety of physical 
phenomena: negative differential resistance Jll, Bloch oscillations [Q, dynamic localization 
0, electric field domain formation [Q, and photon-assisted transport are but a few 
recently studied topics. In the case of weakly coupled multiple quantum well structures the 
perpendicular transport is dominated by resonant transitions between the energy levels in 
neighboring wells. In real systems these resonances are never sharp: they are broadened 
by various scattering mechanisms due to, e.g., impurities, interface roughness, phonons, or 
interactions between the particles themselves (see, e.g., Ref. [0]). For samples with low 
doping impurity scattering may lead to the appearance of impurity bands These ad- 
ditional narrow bands cannot be described by simple Lorentzian spectral functions thus 
necessitating a more detailed microscopic model. In the present paper we develop a fully 
microscopic model for this situation. Our calculations yield a strong temperature depen- 
dence of the current-field relation associated with the occurrence of two different peaks at 
low fields. The implications of these features are compared to both previously published 
and new experimental results. 

For weakly coupled quantum wells the current from the z/*'' level in well n to the /i*'^ level 
in well n + 1 is given by |^ 

dE 



^ J-oo Ztxii 



xA'^^it E) [np^E - fin) - npiE + eFd - fi^+i)] ■ (1) 

Here e is the electron charge, /i„ is the electro-chemical potential in well n, np{£) = 
1/[1 + exp{£ /ksTe)], and Tg is the electron temperature. Fd is the voltage drop per period 
d. The transition matrix elements Hi'" between adjacent wells are calculated using the 
Wannier functions obtained from a Kronig- Penney model |]IU| . Here we restrict ourselves to 
resonant transitions where the tunneling process conserves the wave vector k within the two- 
dimensional electron-gas; nonresonant transitions (see Ref. [|l^) are of minor significance for 
the sample considered here. The microscopical description of the impurity scattering enters 
through the spectral functions A'^{k,E) which are evaluated in a self-consistent single-site 
approximation (shown diagrammatically in the lower inset of Fig. 0) like in Ref. The 
homogeneous doping is modeled by 8 equally spaced 5-doping layers per period. The screen- 
ing of the impurity interaction is accounted for within the random phase approximation. 
We restrict ourselves to the two lowest levels. For the upper level we include a contribution 
to the imaginary part of the self-energy due to optical phonon transitions to the ground 
level, calculated as in |T^. We do not consider the spin-resolved electron-electron interac- 



tion leading to the splitting of the impurity bands for a single impurity location, (the Mott 
transition, see e.g. [|1^]) because the uniform doping will lead to a smearing out of these 
features. The technical details of our calculation will be given elsewhere; here we point out 
that the theoretical framework has recently [jll| been successfully applied to high-doping 



samples, where impurity bands play no role. 

We next describe our numerical results. The upper inset of Fig. |1] shows the spectral 
functions for two different energies below [E = —3 meV) and above {E = 2 meV) the 
bottom of the (free electron) band edge. For E = 2 meV the spectral function exhibits a 
characteristic Lorentzian-type behavior with a peak around E^ = h^k'^ /2m ^ 2 meV. This 
is the generic behavior of a quasiparticle with a finite lifetime due to scattering. In contrast 
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to this behavior the spectral function is quite flat for E = —3 meV, which is a signature of 
an impurity band p. The density of states p{E) = 2/{2tt) J2k ^{k, E) for the lowest level is 
plotted in Fig. |1| (full line). The impurity bands from different locations of impurities overlap 
and manifest themselves in a low density of states in the range -Emin = —5.8 meV< E < 
below the band edge. Importantly, for these low-doped samples the Fermi level is at Ep = —2 
meV, i.e., within the impurity band. As we shall see below, this property has far-reaching 
consequences. 

We calculate the currents In^n+i (eFd) for different electron temperatures by summing 
up the different contributions from Eq. (|Tp. Characteristic results are shown in Fig. ^ for 
an electron density 6 x 10^/cm^ per well provided by the doping. Here we focus in the 
low field region where only tunneling between the lowest levels in each well takes place. 
For Te = 4 K we find a maximum at eFhighrf ~ 5 meV. This is due to tunneling from the 
impurity band to the free states. The maximum occurs at the energy where the bottom of 
the impurity band in one well is aligned with the band edge of the free electron states in 
the neighboring well, i.e., eFhighC? ~ l-f'mini- An increasing temperature leads to a transfer 
of electrons from the impurity band to the free-electron states and consequently the current 
at eFhighci decreases with increasing Tg. The density of states in the impurity band is much 
lower than in the free-electron states, and hence the majority of the electrons will be in the 
free-electron states for ksTe ^ |-Emin|- For the free-electron states the overlap between the 
spectral functions (note that the inset of Fig. |l| is plotted in log-scale) from well n and ra + l is 
very large at small fields. This yields a strong enhancement in the low field conductivity. As 
the overlap of the spectral functions becomes small for voltage differences larger than their 
broadening, the current decreases for eFd > eExo^d 0.5 meV. This transition from a single- 
peaked to a double-peaked current-field relation with increasing temperature is essential in 
understanding the transport in low-doped samples. 

We can now resolve a recent puzzle of Ref . fl^ : there it was found that peaks observed 



in the current-voltage characteristics (IV) under terahertz irradiation can be interpreted 
consistently as photon- replicas of a maximum in an "instantaneous" IV located at ?7 ~ 20 
mV. Yet, a direct measurement of the IV without irradiation showed no maximum at this 
bias. A qualitative explanation, based on impurity bands, is as follows. For low electron 
temperatures and without irradiation the maximum at eFhighC? dominates the domain for- 
mation which sets in at voltages exceeding U NeF^iighd where = 10 is the number 
of wells. If the THz radiation is present the electrons are excited from the impurity band 
into the free electron states corresponding to a larger effective electron temperature. Thus, 
the maximum aX U = NeFio-^d is dominant, and the photon replicas corresponding to this 
feature are seen experimentally. 

While the above consideration explains qualitatively the findings of Ref. [Q, it needs 
to be refined to agree quantitatively with the experiment: the experimental photon- replica 
suggest eF^^d ^ 2 meV which is four times the value from Fig. 0. The domain formation 
sets in a.t U = 0.1 V yielding eF^^^^^d ~ 10 meV (i.e., twice larger than eF^ighd of Fig. |[), 
which seems to require |-Emin| ~ 10 meV. This value follows from our theory as well, as we 
now argue. The screening of the impurity potential was treated using a free electron density- 



of-states po = m/nh'^ yielding the polarizability n(fc) = po 1 — 6{k — 2kF)yl — A{kp/kY 
Now n(0) is related to the actual density of states at Ep ~ —2 meV which is significantly 
lower than pQ. Calculations within the Born- approximation show that the fc-dependence of 
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the polarizability becomes weaker and that n(0) decreases with increasing scattering |T5|,|T6[. 
In order to accommodate these trends we make the replacement Il{k) — >■ Il*{k) = 0.07po, 
given by the calculated density of states of Fig. and use this value in further calculations. 
The smaller polarizability weakens the screening of the impurity potential and therefore 
increases the binding energy of the impurities |p. The resulting density of states is shown 
in Fig. by a dashed line. The onset of the impurity band is now -Emin = —9.1 meV and the 
current-field characteristics look similar to Fig. g but with eFio^d = 2 meV, eFhighrf = 8.2 
meV in good agreement with the values deduced from the experiment. 

As a further test of the theory, we next consider the temperature dependence of the 
zero-bias conductance cr. Eq. (|1|) yields cr{T) oc / dEJ2k^{k,E)^dnp{E)/dE, which implies 
ct(T) oc 1/T for ksT ^ Ep — -Emin and cr(T) — > const, for ksT <ti Ep — -Emm? provided 
that Efc -4(fc, Ef has no strong energy dependence. This is the case if impurity scattering 
is treated within the self-consistent Born approximation where no impurity bands form. 
Then a is monotonously decreasing in T as shown in Fig. |^. However, a different scenario 
emerges if the electrons occupy impurity bands for low temperatures. Then a is strongly 
suppressed due to the small values of A{k, E) foi E < 0, see Fig. |1|. As temperature is 
increased, more electrons are excited to the free electron states, and a increases with T until 
the impurity bands are almost empty at ksT ~ |-Emin|- This physical picture is confirmed by 
our measurements, shown in Fig. ^ together with our full calculation. At low temperatures 
the agreement is quantitative, while at intermediate T the theory overestimates a; this is 
most likely due to additional scattering processes not included in our calculation, or by the 
presence of a contact resistance which may limit the experimental conductance. 

In order to calculate the full current-voltage characteristic we must consider electron 
heating, i.e., estimate the parameter Tg in Eq. (|l]). With average voltage drop of 8 mV per 
well we find experimentally / = 0.57 fiA without irradiation. This gives a power dissipation 
of P = 9 pW per electron and a temperature of Te ~ 47 K using the standard energy loss by 
optical phonons |T^, hui^o/Ti^o exp{—huJi^o/kTe) with tlo ~ 80 fs. An increase of the energy 
loss with decreasing electron density has been predicted due to plasmon-phonon coupling 
||T7| . Recently, P = 9 pW at Te ~ 40 K for a sample with an electron density 6 x 10^°/cm^ 
per period has been reported ||T8|]. As our density is much smaller an even lower electron 
temperature is possible. Experimentally, we find that the current-voltage characteristics 
exhibits a plateau for 0.04 V< U <0.1 V with I = 0.57 fiA (see Fig. The current of 
this plateau remains almost constant for lattice temperatures between 2 K and 30 K, but 
drops for higher temperatures. This can be understood if the electron heating prevents the 
electron temperature from cooling below 30 K. All these facts indicate that Tg ^ 30 K is a 
consistent choice for medium voltages. For vanishing bias, on the other hand, the electron 
temperature must approach the lattice temperature (4K). 

The IV including domain formation [|l9l is calculated as in Ref. [|lTl and is shown in 
Fig. ^. We find excellent agreement with the data for the maximum current for the first 
peak [Ifl^ = 0.8/iA; /f exp = 0.58yuA) and the second peak at higher voltages mih" — 11a*A; 
J^gxp = 13.6/iA, not shown in the figure). The voltage where domain formation sets in is 
in good agreement with the experimental value, but the domain branches occur at higher 
voltages in the experiment. This can be explained by the fact that once the high field domain 
is formed it extends into the low doped receiving contact region where a voltage drop of 0.1 
V may occur. 
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Finally, we consider the sample in a strong terahertz field from a free-electron laser with 
frequency v and field strength F^c- Following Refs. 1^,14-2^ we use the relation 



+n+l 



^ [J/(a)]2/„^„+i(eF,,rf + Ihv) , (2) 
where a = eFacd/ {hv) and Ji is the ordinary Bessel function of order /. This formula was 



originally derived for tunneling through a single barrier |21] and holds within a miniband 
model for superlattices as well [^]. However, modifications due to photon side-bands from a 
single quantum well have been suggested. According to the arguments presented above, 
we model the nonequilibrium distribution function by assuming an electron temperature of 
35 K. Then only 48 % of the electrons are occupying the states below E = Quantitative 
agreement between theory (Fig. ^) and experiment (Fig. ^d) is found for hv = 6.3 meV 



(1.5 THz) for different strengths of the laser field |]2^. We find a direct tunneling peak at 
UdiT = NFiowd ~ 20 mV and photon replicas at f/ ~ Udn + Nhu/e and U ~ t/dir + INhuje. 
For low bias and high intensities there is a region of absolute negative conductance (NC) 
1^ which is studied for different frequencies in the following. In the experiment at each 
frequency the laser power was chosen to give maximum NC near zero bias. As can be seen 
from Fig. ^(d) a minimum in current occurs at f/ ^ — f^dir + Nhv/e. This is just the photon 
replica of the direct tunneling peak on the negative bias side which dominates if the 
direct tunneling channel is suppressed close to the zero of Jq^o) in Eq. (|^), i.e., a ~ 2.4, 
as used in the calculation (Fig. ^). For hv = 5.3 meV a smaller value of a (thin line) 
agrees better with the experimental data. This may be related to charge accumulation 
effects inside the sample weakening the NC, so that maximum NC is observed at a laser 
field corresponding to a lower value of a. Both the theoretical and experimental result show 
that absolute NC is persistent in a wide range of frequencies but becomes less pronounced 
with decreasing photon energy. 

In conclusion, we have presented a self-contained theory of transport in low doped and 
weakly coupled quantum wells. A detailed comparison with several different experiments 
shows that: (i) it is necessary to consider impurity bands; (ii) excellent agreement between 
theory and experiment is achieved for static and irradiated IV as well as zero-bias con- 
ductance using reduced screening, and (iii) the simple relation (0) provides a quantitative 
description of photon-assisted transport in multiple quantum wells. 
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FIG. 1. Calculated density of states in units of the 2D free carrier density po using screening 
due to a free electron gas (full line) and reduced screening (dashed line). The upper inset shows 
the spectral function A{Ek,E) for different energies E. The lower inset depicts diagrammatically 
the self-consistent single-site approximation. 
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FIG. 2. Current-field relation for different electron temperatures. 
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FIG. 3. Comparison of experimental and theoretical conductance at zero bias versus temper- 
ature. 
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FIG. 5. Current-voltage characteristics under irradiation, a) Theoretical results for hv = 6.3 
meV and different field strength eFacd = ahu of the irradiation, b) Experimental results for 
hf = 6.3 meV and different laser intensities increasing from the top to the bottom. The actual 
values Fac inside the sample are not accessible, c) Theoretical results for a = 2.4 and different 
photon energies. The thin line depicts hv = 5.3 meV and a = 2.1. d) Experimental results for 
different photon energies. The laser intensity was chosen to give maximum negative conductance. 
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